Indications on neutrino oscillations parameters 
from initial K2K and current SK data 



G.L. Fogli, E. Lisi, and A. Marrone 

Dipartimento di Fisica and Sezione INFN di Bari 
Via Amendola 173, 70126 Bari, Italy 



Abstract 



We briefly discuss the impact of initial data from the KEK-to-Kamioka (K2K) 
neutrino experiment on the z/^ — > Vr oscillation parameters {m? ,ta.v? ijj) cur- 
rently indicated by the Super-Kamiokande (SK) atmospheric neutrino exper- 
iment. After showing the very good agreement between K2K and SK, we 
combine the two separate pieces of information. We find that the 99% C.L. 
range for m? allowed by SK only, m? E [1.3, 5.6] x 10~^ eV^, is reduced to 
[1.5, 4.8] X 10^^ eV^ by including K2K data. By halving the uncertainties of 
the K2K total rate (with central value unchanged), the m? range would be 
ulteriorly reduced to [1.8, 4.0] x 10~^ eV^. Such information appears to be al- 
ready useful in planning (very) long baseline neutrino oscillation experiments. 
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I. INTRODUCTION 



The KEK-to-Kamioka (K2K) long baseline neutrino experiment |^ is designed to explore, 
with laboratory neutrinos, the mechanism of t'^ disappearance indicated by the Super- 
Kamiokande (SK) 0, MACRO 0, and Soudan 2 |^ atmospheric neutrino experiments. 
Such mechanism appears to be dominated by z/^i Vr oscillations f^-QI with mass-mixing 
parameters (m^, tan^ close to (3 x 10~^ eV^, 1). The pathlength (L = 250 km) and typ- 
ical energies {E ~ few GeV) of K2K i/^'s are well suited to study disappearance effects 
for such parameter values |l|]. 

The recent data released by the K2K experiment [P|-pT| (44 observed neutrino events 
vs 63.9 expected) are already inconsistent with the no-oscillation hypothesis at about 97% 
CL. [0-0, and suggest disappearance of the muon neutrino flavor in the path from KEK 
to Kamioka. The K2K collaboration is being understandably rather conservative on the 
oscillation explanation of the data, pending more accurate estimates of the K2K neutrino 
spectrum and interaction uncertainties by means of the near detector |T^. However, it is 
tempting to use just the minimum amount of spectrum-integrated data (i.e., the total event 
rate, for which the error estimate appears to be relatively stable), in order to study the 



compatibility of K2K and SK data [|14 



In this paper we attempt to investigate such a compatibility issue, by comparing predic- 
tions and data in both K2K and SK, within the simplest (two-family) — v^. oscillation 
framework. We show that the K2K and the SK data are consistent, and that their combi- 
nation starts to be useful in constraining further the oscillation parameters. In particular, 
the 99% CL. range for m? allowed by SK only, im? G [1.3, 5.6] x 10~^ eV^, is reduced to 
[1.5, 4.8] X 10~^ eV^ by including K2K data. Finally we show that, by halving the K2K 
total rate error, the rn? range would be ulteriorly reduced to [1.8, 4.0] x 10~^ eV^. 



II. ANALYSIS OF K2K 

Analyses of K2K data have to consider the energy spectrum of parent neutrinos at the 
far detector S{E), i.e., the spectrum of neutrino that have produced a detected interaction 
at the Kamioka site. Such spectrum is, in general, given by the product of the unoscillated 
neutrino energy spectrum $ at the far (SK) detector, times the interaction cross section 
a, times the detection efficiency e (all in differential form), integrated over all the final 
state parameters (simbolically, S = /$ ■ cr ■ e). The spectrum S{E) is important since, 
in the presence of neutrino mass and mixing, the oscillated event rate can be expressed as 
r = / dE S{E) ■ P^^{E) [to be compared with the unoscillated rate, tq = / dE S{E)], where 
Pfj,f^ is the z/^ survival probability. 

At present, there is enough public information on $ |^11|, but not yet on the product a-e 



relevant for the far detector specifications. Therefore, a direct and accurate reconstruction 
of the spectrum S{E) and of its uncertainties is not possible outside the K2K collaboration, 
at least at present. However, an indirect, approximate reconstruction can be made, by 
considering that: (a) the neutrino spectrum is practically zero above 5 GeV [P]; and (b) 
according to the K2K MonteCarlo (MC) simulation, for the current exposure (3.85 x 10^^ 
p.o.t.) one has tq = 63.9 events, while r = 41.5, 27.4, and 23.1 events for = 3, 5, and 
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7 X 10~^ eV^, respectively, assuming maximal mixing p|-pl|. We have then introduced an 
empirical functional form for the spectrum, S{E) oc — x)^ with x = E/{5 GeV), which 
is found to reproduce well the above MC results for (a, /?) = (1.50, 3.34). In the following 
analysis we use such interacted neutrino spectrum, since it gives — by construction — results 
close to the K2K MC simulation, after integration over P^^{E). A more accurate (and less 
tentative) parametrization of S{E) shall be possible after the experimental specifications 
and selection cuts (relevant to determine a ■ e and its uncertainties) will be measured and 
described in detail. 

Concerning the uncertainties, the current systematic error on the total rate is estimated 
to be 10% by the K2K Collaboration p|^TT|. Spectral bins may be affected by larger (and cor- 
related) errors |]9HTI|, which, however, do not affect our spectrum-averaged analysis. There- 
fore, we attach a 10% fractional error to the oscillated number of events r = r(m^, tan^ V^), 
while the statistical (Poisson) error cistat = "\/44 is attached to the number of observed 
events, rgxp = 44. The statistics for K2K is then simply given by (r — rexp)^/^"^^^, where 
^tot = ^stat + (0.1 r)^. With such definition, the no oscillation hypothesis (r = tq) is disfa- 
vored at 97% C.L. (x^ = 4.7 for A^df = 1), consistently with the current claim of the K2K 
Collaboration IHg. 

Figure 1 shows the number r of K2K events for the current exposure (3.85 x 10^^ p.o.t.), 
as a function of m?, for maximal mixing {ian^ ip = 1). The solid curve, corresponding to 
our calculation of r{iv?), starts from 63.9 events in the no-oscillation limit {m? 0), and 
tends to the asymptotic value of 63.9/2 in the fast oscillation limit {m? — * oo), after a few 
"wiggles" associated to the first oscillation cycles. The curve is very close to the four K2K 
MC points used to benchmark our calculation. The horizontal gray band represents the 
ilcTtot interval for the total number of observed events in K2K. The allowed band disfavors 
no oscillations at 2.2o", and the first deep oscillation minimum at 2.7cr. It is instead in 
very good agreement with the oscillated predictions for im? ~ few X 10^^ eV^, the best fit 
being reached at rn? = 2.7 x 10~^ eV^. Such a value is very close to the one independently 
quoted by the SK Collaboration as the current best fit to their atmospheric neutrino data 



(m 



^ — 2.5 X 10 ^ eV^) [|T5l, under the same z/^ Ur oscillation hypothesis.^ Therefore, it 



makes sense to combine K2K and SK data, as we do in the following section. 



III. K2K AND SK: COMBINATION AND AND IMPLICATIONS 

Given the very good agreement between the K2K and SK allowed ranges for at 
maximal mixing, it is interesting to study their compatibility for unconstrained 2i> mixing. 
Concerning the SK atmospheric 2z/ analysis, we make use of the results recently reported 
by us in p (for 79.5 kTy SK data), specialized to the simplest scenario of pure t'^ — > Ur 
oscillations.0 



^It is also very close to the best fit of upgoing muon data in MACRO, m? = 2.4 x 10 ^ eV^ Q. 

^Therefore, in this work we use A^df = 2 to draw iso-Ax^ contours, while we used Adf = 3 for 
the more general 3z/ and 4z/ cases considered in [^. 
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Figure 2 shows the regions allowed at 90% and 99% CL. by SK, K2K, and their combina- 
tion, in the 2v mass-mixing plane (m^, tan^ ip). The upper left panel corresponds to SK atmo- 
spheric data only, for which we find two degenerate best-fit points (stars) at m? = 3 x 10~^ 
eV^ and at octant-symmetric mixing values, {ian^ ip)^^ = 0.76, corresponding to slightly 
nonmaximal oscillation amplitude, sin^ 2ip = 0.98.0 

The upper right panel in Fig. 2 shows the fit to K2K only. The locus of best-fit points 
is a continuous, octant-symmetric curve (not shown) passing through (m^/eV^, tan^ ip) = 
(2.7 X 10~^, 1). The K2K constraints in the mass-mixing plane are weaker than those 
placed by SK, especially on tan^ ip. Therefore, one cannot expect a significant improvement 
on tan^?/' limits from the SK+K2K combination. Concerning m^, the no-oscillation limit 

— > is still allowed at 99% C.L., while values around ~ 8.5 x 10~^ eV^ are excluded 
(for large mixing). Such values correspond to the first (deep) oscillation minimum in Fig. 1, 
which is disfavored by the data at 2.7a, as previously noted. Therefore, in the SK+K2K 
combination, we expect "high" values of to be more disfavored than "low" values. 

The lower left panel in Fig. 2 shows the combination of SK and K2K data. The best-fit 
points are located at the same values of tan^ ip as for SK alone and, in general, the bounds 
on tan^ ip are not significantly modified, as expected. The best-fit value of is only slightly 
lowered (m^ = 2.9 x lO""^ eV^) but, most importantly, the 90% and 99% CL. ranges of 
are appreciably reduced both from below and (more strongly) from above. Such results show 
that the K2K experiment is already having a nonnegligible impact in the determination of 
the neutrino squared mass difference relevant for the z/^ Vr channel and for the leading 
oscillations in (very) long baseline experiments. 

The lower right panel in Fig. 2 shows a prospective SK+K2K combination, with the 
same SK data and K2K data but with K2K total uncertainty hypothetically reduced by a 
factor of two. The iv? range is significantly narrowed, while the tan^ ip range is still basically 
unchanged. In summary, the various panels of Fig. 2 demonstrate that K2K data are relevant 
for the determination of m^, rather than of i&v? ip. Therefore, it makes sense to discuss in 
more detail the results of Fig. 2 in terms of only, with tan^ ip unconstrained. 

Figure 3 shows the results of such an exercise, in terms of the function x^(m^). The 
90% and 99% CL. ranges for (A^df = 2) are explicitly shown for the SK and SK+K2K 
fits. Numerically, we find that the 99% CL. range for allowed by SK only, G 
[1.3, 5.6] X 10~^ eV^, is reduced to [1.5, 4.8] x 10~^ eV^ by including K2K data. By halving the 
uncertainties of the K2K total rate, the m? range would be ulteriorly reduced to [1.8, 4.0] x 
10-3 eV2 (99% CL.). 

The reduction of the minimum values of th? in Fig. 3 (at a given CL.) is relevant for 
future (very) long baseline experiments. For instance, in the OPERA experiment |]T^ down- 
stream the CERN-to-Gran Sasso neutrino beam, the r appearance rate is approximately 
proportional to (m^)^ l]!^; therefore, the increase of w^mm froni 1.3 x 10"^ eV^ (SK) to 
1.5 X 10"^ eV^ (SK+K2K) at 99% CL. implies an increase of the minimum expected r 
event rate by a factor ~ 1.3 (at the same CL.); the increase would be as high as a factor 



^However, such small deviations from maximal mixing at best fit (also found in [|T6|) are not 
statistically significant at present. 



4 



~ 2 for the case with halved K2K errors (m^jj^ = 1.8 x 10~^ eV^). The reduction of the 
allowed range is also important to refine the energy or baseline optimization in future neu- 
trino factory experiments, where the current uncertainty plays an important role (see, 
e.g., |jl8| and references therein). Finally, the increase of m^;^ strengthens the accuracy of 
approximations based on the hierarchy of squared mass differences, which have long been 



used in global solar+terrestrial neutrino analyses (see, e.g., [g,|Tg,|I^). 

A final remark is in order. Our current analysis is based on initial K2K data and on an 
approximate description of the K2K detector specification. Therefore, it cannot be a substi- 
tute of the (joint) official oscillation analysis that will be performed by the K2K (SK+K2K) 
collaboration(s). However, we think that our results, although necessarily approximate, are 
sufficiently indicative of the K2K potential in improving our current knowledge of z/^ Vr 
oscillations. 



IV. SUMMARY 



We have shown that, within the simplest (two-family) i/^ — oscillation scenario, the 
initial evidence for a flux suppression in K2K is perfectly consistent with the atmospheric 

flux anomaly. The range of the neutrino squared mass difference indicated by the SK 
atmospheric v experiment is reduced by including the K2K data. As a consequence, K2K 
starts to be important in narrowing the range of predictions for future (very) long baseline 
experiment. 
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FIG. 1. Number of events in K2K for the current exposure (3.85 x 10^^ p.o.t.), as a function 
of m? (at taia^ ip = 1). Our calculation (solid curve) is benchmarked by the K2K MC simulation 
(square markers). The horizontal gray band represents the current K2K data within one standard 
deviation. 
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SK + K2K, 2v oscillations 




FIG. 2. Separate two-flavor oscillation analyses of SK and K2K, together with their combination 
(with eventual halving of K2K errors). See the text for details. 
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Bounds on for unconstrained 2v mixing 




FIG. 3. Bounds on at 90% and 99% C.L. (ATdf = 2) from the two-flavor analysis of the 
SK and K2K data. 



9 



